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Abstract 

The calculation of the next-to-leading order QCD corrections to the pro- 
duction of a neutral Higgs-boson pair is briefly summarized. The dominant 
production mechanism is gluon fusion via top-quark loops, so that the QCD 
corrections can be considered, to good approximation, in the limit of a heavy 
top-quark mass. 
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1. Introduction 

The search for Higgs bosons and the investigation of their properties will 
be one of the most important tasks in future high-energy collider experi- 
ments. Direct information on the Higgs potential, which reveals the details 
of the mass-generation mechanism in spontaneously broken gauge theories, 
can only be obtained by measuring the self-interactions of the found Higgs 
boson(s), i.e. by searching for multiple-Higgs production processes. In this 
article we briefly summarize the results of Ref. Q , in which we discussed the 
QCD corrections to the production of neutral Higgs-boson pairs in pp colli- 
sions at next-to-leading order (NLO) in the framework of both the Standard 
Model (SM) and its minimal supersymmetric extension (MSSM). 

The most important production mechanism of Higgs-boson pairs in pp 
collisions is gluon fusion through a heavy-quark loop, i.e. a NLO prediction 
requires a two-loop calculation. The calculational standard for multiloop 
diagrams is, however, not high enough yet to allow for an exact analytical 
evaluation of massive two-loop vertex and box diagrams. Fortunately, the 
structure of the relevant diagrams admits a feasible approximation that is 
also phenomenologically reasonable. It consists in an asymptotic expansion 
in a heavy top-quark mass, rrit — > oo, and the neglect of the 6-quark mass, 
mb = 0. Therefore, we have calculated the relative NLO corrections in this 
limit, but taking into account the known leading-order (LO) cross section 
[0 exactly. The analogous procedure for single-Higgs production at LHC 
energies reproduces the exactly known two-loop result within 5% for Higgs- 
boson masses below the 2mt threshold of the top-quark loops; even for 
Higgs-boson masses up to 1 TeV the deviation does not exceed 10%. This 
impressive reliability is due to the dominance of universal soft-gluon effects 
, which do not resolve the detailed structure of the hard-scattering process 
and are fully included in the approximation. Therefore, we expect that our 
results for Higgs-pair production are valid at the 10% level, at least for 
invariant masses of the Higgs-boson pair below 2mt. 

For the production of Higgs-boson pairs in the MSSM a few more re- 
marks are in order. The neglect of forces us to choose tg/3 small, since 
the 6-quark Yukawa coupling is enhanced, and thus non-negligible, for large 
tg/3. Moreover, we consistently neglect squark-loop contributions, i.e. we 
assume squark masses to be large. The production channels for a scalar 
and a pseudoscalar Higgs-boson pair, hA and HA, represent an exceptional 
case, because they receive also non-negligible contributions from the Drell- 
Yan-like mechanism via qq fusion, which are, however, well known at LO 
and NLO. More details can be found in Ref. 
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2. Calculational framework 

The LO predictions for the gluon- fusion processes gg — > (pi(j)2, where 
<j)i(j)2 stands for any neutral Higgs-boson pair, via quark loops are known, 
including the exact dependence on the heavy-quark masses The follow- 
ing one-loop diagrams contribute: 




where the Z-boson exchange in the s-channel is only relevant for hA and 
HA production. At NLO, we get virtual corrections formed by two-loop 
diagrams, which contain a closed quark loop and one, two or three internal 
gluons. Moreover, there are real one-loop corrections induced by the parton 
processes gg — > (j)i(j)29, 9Q — > 4'i4'2Q, and qq — > (j)i(j)29- After summing vir- 
tual and real corrections, IR divergences drop out; the remaining collinear 
singularities are absorbed in the parton distributions when calculating the 
pp cross section. 

For the heavy-mass expansion of the diagrams containing the top quark 
we employ the general algorithm of Ref . . The actual expansion is carried 
out in the integrand of a diagram, i.e. before the momentum-space inte- 
gration, and relies on dimensional regularization. The general algorithm is 
easily summarized. Given any Feynman graph F with some large internal 
masses Mj, and denoting the corresponding amplitude and integrand by Fy 
and Ir, respectively, the large-mass expansion reads 

Ft = I{Uid^qi)Ir KJlid^qi) Ir/.Tj^im^, (1) 

where qi are the integration momenta. The sum on the r.h.s. runs over all 
subgraphs 7 of F that contain all propagators with the heavy masses Mj 
and that are irreducible with respect to those lines of 7 that carry light 
masses mj. The integrand of the subgraph 7 is denoted by I^. The reduced 
graph F/7 results from F upon shrinking 7 to a point, and the integrand 
/r/7 is defined such that Ir = I^yl-p/^. The symbol Tpi represents an 
operator that replaces the integrand /-y by its Taylor series in the expansion 
parameters pj and mj, where pj are the external momenta of the subgraph 
7. Therefore, Eq. (|^) expresses the original integral i*r by a sum over simpler 
integrals with a finite number of terms in each order M". 

As an illustrating example, we inspect the following box diagram. 



There are two subgraphs 7 that are relevant in the expansion (||) of this 
graph. The first one, 7ioop) consists of all internal lines, the second one, 7top, 
consists of the top-quark lines only. The reduced graphs look as follows. 




r/7ioop: m r/7, 




7top- 



g r \ (p2 g 

In the case of 7ioop the Taylor-expansion operator replaces each propagator 
P{q — p, m) with m = 0, mj by 

P{q -p,m) = [{q - pf - m^]"! = E/=o(9^ " m'^y^-\2qp - p^)\ (2) 

where g is a linear combination of the two integration momenta qi and 52, 
and p consists of external momenta pi . This replacement leads to two- loop 
vacuum integrals of the form 

'"^ V {qlY^ {ql - ml)- [{q, + q,Y - mfj-a ' ^"^^ 

the calculation of which is straightforward in n dimensions. For 7top the 
Taylor expansion concerns only the top-quark propagators, and the replace- 
ment (H) applies for m = rrit and q = q2, where q2 is the loop momentum 
running through the top-quark loop. Note that p now includes all external 
momenta of the process as well as the loop momentum qi running through 
the internal gluon lines. Thus, the integration over (72 yields simple one-loop 
vacuum integrals 

while the q\ integration leads to one-loop integrals of the form 



(5) 



containing massless propagators and non-vanishing external momenta pi 
and p2- 
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It is interesting to realize the close connection of the diagrammatic ex- 
pansion algorithm described above and the independent approach using a 
low-energy effective theory. The effective Lagrangian results from integrat- 
ing out the heavy top quark from the full theory and quantifies the local 
interactions between light Higgs bosons and gluons at low energies induced 
by heavy-quark loops. The effective interactions that are relevant in our 
case are known in LO as well as in NLO and are summarized in Ref. Q. 
The effective couplings can be derived from the gluon self-energy upon dif- 
ferentiation, except for the ones that involve an odd number of pseudoscalar 
Higgs bosons. The latter are related to the ABJ anomaly. We have also 
calculated the NLO corrections to gg — > (j)i4'2 in the effective-Lagrangian 
approach. There is a one-to-one correspondence between the sum of all 
topologically equal reduced graphs r/7 in the diagrammatical expansion 
and the corresponding effective diagram in which 7 is identified with the 
effective coupling. More precisely, the effective couplings enter the NLO 
corrections in two different ways. Firstly, the NLO parts of the effective 
interactions contribute in tree- level diagrams like r/7ioop- Secondly, the LO 
parts of the effective couplings contribute in one-loop diagrams like r/7top, 
which are the only sources of non-local effects. 

Although we set = 0, there are contributions from 6-quark loops in 
the hA and HA production channels. These contributions are associated 
with Z-boson exchange and have to be considered separately. The virtual 
two-loop corrections in the Zgg vertex are fully detemined by gauge invari- 
ance, because only a single form factor contributes. The divergence of the 
Zgg vertex is related by a Ward identity to the coupling of a pseudoscalar 
field to two gluons, implying that the Zgg vertex correction concides with 
the one for single A production, which gets no 6-quark contribution for 
mi, = 0. The 6-quark loops appearing in the real corrections do not possess 
a simple form. However, the dominant effects, the ones that are sensitive 
to soft gluons, are already included by taking into account the Zgg vertex 
only, leading again to the same correction as for the Agg vertex. 



3. Results 

The final analytical results for the NLO corrections to the process pp 
0102 + ^ have a rather compact form in the adopted approximation. In 
this context, the LO cross section (Tlo of the parton process gg — > 0i02 is 
of central importance: 

aLo{S = Q^) = Jdi^^^{\C^FA + CaFaf + \CaGa\'}, (6) 
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where Q is the invariant mass of the Higgs-boson pair, and s, i, u are the 
usual Mandelstam variables in the partonic centre-of-mass (CM) system. 
The functions Ca.d contain coupling factors and s-channel propagators, and 
-Fa,d and Gn represent the form factors corresponding to the two tensors 
that are relevant for spin-0 and spin-2 exchange. The explicit expressions 
of the Cs, -F's, and G's, including the full fermion-mass dependence for 
the latter two, are given in Ref. [^. Here we only give the Cs and the 
asymptotic form of the form factors for the SM: 

_ 3M^ _ 

Fa ^ |, Fn^-|, Gn ^ for ^ oo, m;, = 0. (7) 

Note, however, that the exact form of the form factors is used in the nu- 
merical evaluations. 

All contributions to the pp cross section cjnlO) 

can be written as a convolution of a correction factor, of o"lO) and of the 
parton-parton luminosities dC^^ /dr: 

o-LO = / dr — — o-lo(ts), 
Jro dr 

AfJvirt = / dr — — fJLo (ts) C, 

vr Jro dr 



n ^ / / T '"^^^'^'^ -^P..(^)log — 



2„ , M2 2 



^ log + - (1 - ^)^| , 

where i/s is the total CM energy in the pp system, and y/r^s = mi +7712 is the 
corresponding threshold for the production of two Higgs bosons with masses 
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123.5 131.6 143.3 157.9 174.4 192.0 210.3 Mjj [GeV] 
67.1 77.7 85.0 89.6 92.4 94.3 95.5 MJGeV] 



2 - 



K(pp (j)i(j)2+X) 
Vs=14 TeV 



m, = 175 GeV 



hA 



HA 




AA 



qq HA 
qq hA 

J 1 1 1 L 



80 



100 



120 



140 
M, [GeV] 



160 



180 



200 



Fig. 1. K factors of the QCD- corrected gluon-fusion and Drell-Yan like cross 
sections a(pp <pi(f>2 + X) at the LHC with CM energy ^/s = 14 TeV. 



mi and m2- The renormalization and factorization scales are denoted by 
and M, respectively, and Pij{z) are the Altarelli-Parisi splitting functions. 
While the form of the real corrections Aagg^gg^gq is universal, the virtual 



correction factors 

C = vr^ + ci + log ^ (10) 

/ di |c2 Cn(CAFA + CuFu) + ^^^P (f + f ) CgG: 



+ 3f?e 



Jdi{\CAFA + CnFn|2 + \CnGn\^} 



do not coincide for all Higgs-boson pairs. The coefficients Cj for the individ- 
ual final-state Higgs bosons 0i</>2 are given by 

4ii<j)2 = hh, hH, HH : ci = ^, C2 = |, C3 = |, ^4 = |, 
(j)i(j)2 = hA, HA : ci = 6, C2 = |, 03 = -|, C4 = |, 

^l4>2 = AA: ci = ^, C2 = -1, C3 = 1, C4 = 1. (11) 



A detailed numerical study of the QCD corrections described above is 
presented in Ref. [|| for the SM and the MSSM. Figure |l| summarizes our 
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results on the K factors for Higgs-pair production in the MSSM as func- 
tions of the pseudoscalar mass M^, demonstrating their phenomenological 
importance at the LHC. The K factor for the SM cross section smoothly 
varies between 1.9 and 2.0 for a Higgs-boson mass within 80-200 GeV. 
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